(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property 
Organization 

International Bureau 

(43) International Publication Date 
29 July 2004 (29.07.2004) 




PCT 



(10) International Publication Number 

WO 2004/062464 A2 



(51) International Patent Classification 7 : 



A61B 



(21) International Application Number: 

PCT/KR2003/002825 



(22) International Filing Date: 

24 December 2003 (24.12.2003) 



(25) Filing Language: 

(26) Publication Language: 



English 
English 



(30) Priority Data: 

10-2003-0001622 10 January 2003 (10.01.2003) KR 
10-2003-0041569 25 June 2003 (25.06.2003) KR 

(71) Applicant (for all designated States except US): KO- 
HWANG FOUNDATION, KOHWANG BOARD OF 
TRUSTEE [KR/KR]; Hoegi-dong 1, Dongdaemun-gu, 
130-701 Seoul (KR). 



(71) Applicants and 

(72) Inventors: WOO, Eung Je [KR/KR]; Samick Villa 
106-104, Gumi-dong 121, Bundang-gu, Seongnam-si, 
463-743 Gyeonggi-do (KR). KWON, Ohin [KR/KR] ; 
Hanshin-Chunggu APT., 112-1502, Mok 6-dong, 
Yangcheon-gu, 158-751 Seoul (KR). SEO, Jin Keun 
[KR/KR]; Samhogarden Mansion Ra-308„ Banpo-dong 
30-18, Seocho-gu, 137-040 Seoul (KR). 

(74) Agents: BAHNG, Hae Cheol et al.; Kims International 
Patent & Law Office, 15th Floor Yo Sam Building, 648-23, 
Yeoksam-dong, Kangnam-gu, Seoul 135-080 (KR). 

(81) Designated States (national): AE, AG, AL, AM, AT, AU, 
AZ, BA, BB, BG, BR, BW, BY, BZ, CA, CH, CN, CO, 
CR, CU, CZ, DE, DK, DM, DZ, EC, EE, EG, ES, FI, GB, 
GD, GE, GH, GM, HR, HU, ID, IL, IN, IS, JP, KE, KG, 
KP, KZ, LC, LK, LR, LS, LT, LU, LV, MA, MD, MG, MK, 
MN, MW, MX, MZ, NI, NO, NZ, OM, PG, PH, PL, PT, 
RO, RU, SC, SD, SE, SG, SK, SL, SY, TJ, TM, TN, TR, 
TT, TZ, UA, UG, US, UZ, VC, VN, YU, ZA, ZM, ZW. 

(84) Designated States (regional): ARIPO patent (BW, GH, 
GM, KE, LS, MW, MZ, SD, SL, SZ, TZ, UG, ZM, ZW), 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 

[ Continued on next page ] 



(54) Title: SYSTEM AND METHOD FOR THREE-DIMENSIONAL VISUALIZATION OF CONDUCTIVITY AND CURRENT 
DENSITY DISTRIBUTION IN ELECTRICALLY CONDUCTING OBJECT 




< 
^© 

>sO 

(57) Abstract: System for visualizing conductivity and current density distributions including a plurality of current injecting devices 
for injecting currents into a measuring object, an MRI scanner for measuring one directional component of a magnetic flux density 
due to each of the currents injected into a measuring object, an operating part for controlling the current injecting devices so as 
to inject currents of different directions into the measuring object, and calculating a conductivity distribution and a current density 
distribution inside of the measuring object by using the one directional component of a magnetic flux density, and displaying means 
for visualizing the conductivity and current density distributions calculated by the operating part, thereby permitting to visualize the 
conductivity and the current density of the measuring object more accurately. 



WO 2004/062464 A2 II I II I « 



European patent (AT, BE, BG, CH, CY, CZ, DE, DK, EE, For two-letter codes and other abbreviations, refer to the "Guid- 
es, FT, FR, GB, GR, HU, IE, IT, LU, MC, NL, PT, RO, SE, ance Notes on Codes and Abbreviations" appearing at the begin- 
SI, SK, TR), OAPI patent (BF, BJ, CF, CG, CI, CM, GA, ning of each regular issue of the PCT Gazette. 
GN, GQ, GW, ML, MR, NE, SN, TD, TG). 

Published: 

— without international search report and to be republished 
upon receipt of that report 



WO 2004/062464 



1 



PCT/KR2003/002825 



SYSTEM AND METHOD FOR THREE-DIMENSIONAL 
VISUALIZATION OF CONDUCTIVITY AND CURRENT 
DENSITY DISTRIBUTION IN ELECTRICALLY CONDUCTING 

OBJECT 

5 

Technical Field 

The present invention relates to system and method for visualization of 
conductivity (or resistivity) and/or current density distributions of a measuring object, 
10 such as a human body, or a substance. 

Background Art 

In general, X-ray, MRI, or ultrasound is used in visualization of an inside 
structure of the human body or the substance. However, the methods can not 
15 visualize electrical properties of the human body or the substance. 

As an effort for solving the problem, a method for visualizing current density 
distribution of an inside of a measuring object using the MRI is initially suggested by 
a research team of the Toronto University in 1989, and, thereafter, there have been 
active researches related thereto. This is the Current Density Imaging (CDI) 
20 technique. 

In the present CDI, a magnetic flux density B due to an injected current I is 
measured by using MRI technique, and a current density J is calculated by using the 

Ampere's law J = — VxB, for visualizing an inside current density distribution. 
Mo 

However, the CDI has a drawback in that rotation of the measuring object in 
25 an MRI scanner is required for obtaining the magnetic flux density B having three 
components Bx, By, and Bz. This is because the MRI scanner can only measure a z- 
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directional component, a direction the same with a main magnetic field, of the flux 
density, i.e., Bz, when the measuring object is in the MRI scanner. 

That is, since the component the MRI scanner can measure at a time is only 
the Bz component, for obtaining all the required three components of the magnetic 
5 flux density vector, the present CDI technique has a serious drawback of requiring 
the measuring object (the human body, or the substance) to be rotated in the MRI 
scanner. 

Meanwhile, as a known method for visualizing the electrical properties of 
the human body, or the substance, there has been Electrical Impedance Tomography 
10 (EIT) that has been under active research starting from late 1970s. The EIT provides 
an image of resistivity (or conductivity) distribution, an electrical property of the 
measuring object. 

In the EIT mostly taking the human body as the measuring object, many 
electrodes are attached on a surface of the human body, for visualizing the resistivity 
15 distribution of an inside of the human body. 

The visualization of the human body according to the resistivity is made 
possible since tissues of the human body, such as blood, bones, lung, heart, and the 
like have electrical properties different from one another. 

However, due to a fundamental drawback of the EIT, an image of the EIT is 
20 poor, to support only a low resolution. That is, EIT has a fundamental drawback in 
that a current-voltage data measured by EIT is extremely insensitive to variation of 
resistivity of an inside of the human body. Therefore, clinical application of EIT is 
difficult, presently. 

Even if CDI is applied to EIT, because it is required to rotate the measuring 
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object (the human body, or the substance), the technical drawback can not be 
resolved. 

Disclosure of Invention 
5 An object of the present invention, designed to solve the related art problem, 

is to provide system and method for visualization of conductivity and current density 
distribution, in which electrical properties (conductivity and current density) of a 
measuring object can be visualized in more accurate high resolution without rotating 
the measuring object. 

10 The object of the present invention can be achieved by providing a system 

for visualizing conductivity and current density distributions including a plurality of 
current injecting devices for injecting currents into a measuring object, an MRI 
scanner for measuring one directional component of a magnetic flux density due to 
each of the currents injected into a measuring object, an operating part for controlling 

15 the current injecting devices so as to inject currents of different directions into the 
measuring object, and calculating a conductivity distribution and a current density 
distribution inside of the measuring object by using the one directional component of 
a magnetic flux density, and displaying means for visualizing the conductivity and 
current density distributions calculated by the operating part. 

20 The current injecting device includes a current source and plurality of 

recessed electrode assemblies. Each recessed electrode assembly includes an 
electrode, an insulating container with the electrode attached at one side and an 
electrolyte substance inside, and a wire for supplying the current to the electrode. 

The operating part calculates an inside voltage and a surface voltage of the 
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measuring object for a first conductivity, and calculates a second conductivity by 
using the inside voltage and the measured one directional component of the magnetic 
flux density. The operating part multiplies or divides a constant to the second 
conductivity according to a ratio of a calculated surface voltage to a measured 
5 surface voltage. 

The operating part determines that the second conductivity is a true 
conductivity, if an absolute value of a difference of the first conductivity and the 
second conductivity is smaller than a tolerance value, and if the absolute value of the 
difference of the first conductivity and the second conductivity is greater than the 
10 tolerance value, the operating part calculates the inside voltage and the surface 
voltage for the second conductivity, and calculates a third conductivity by using the 
inside voltage for the second conductivity and the measured one directional 
component of the magnetic flux density. This process repeats until the tolerance 
value is achieved. 

15 In other aspect of the present invention, there is provided a method for 

visualizing conductivity and current density distributions including the steps of (a) 
injecting currents of different directions into a measuring object through current 
injecting devices attached to a surface of the measuring object, (b) measuring one 
directional component of a magnetic flux density of due to each of the currents 

20 injected into the measuring object, (c) calculating a conductivity and a current 
density of the inside of the measuring object by using the measured one directional 
component of the magnetic flux density, and (d) visualizing the conductivity and the 
current density. 

The step (a) includes the steps of selecting a pair of the current injecting 
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devices in succession, and injecting the currents into the measuring object in 
succession through the selected pair of the current injecting devices. 

The step (c) includes the steps of calculating the inside voltage and the 
surface voltage of the measuring object for a first conductivity, and calculating a 
5 second conductivity by using the calculated inside voltage and the measured one 
directional components of the magnetic flux density. 

The step (c) further includes the step of multiplying or dividing the second 
conductivity by a constant according to a ratio of the calculated surface voltage and a 
measured surface voltage of the measuring object, and further includes the steps of 

10 calculating an absolute value of a difference of the first conductivity and the second 
conductivity, and comparing the absolute value of the difference to a tolerance value. 

The step (c) further includes the steps of determining the second 
conductivity being a true conductivity, if the absolute value of the difference of the 
first conductivity and the second conductivity is smaller than the tolerance value, and 

15 opposite to this, if the absolute value of the difference of the first conductivity and 
the second conductivity is greater than the tolerance value, the step (c) further 
includes the steps of calculating the inside voltage and the surface voltage of the 
measuring object for the second conductivity, and calculating a third conductivity by 
using the inside voltage and the measured one directional component of the magnetic 

20 flux density for the second conductivity. 

Brief Description of Drawings 

The accompanying drawings, which are included to provide a further 
understanding of the invention, illustrate embodiment(s) of the invention and 
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together with the description serve to explain the principle of the invention. In the 
drawings; 

FIG. 1 illustrates a diagram of a system for visualizing conductivity and 
current density distribution in accordance with a preferred embodiment of the present 
5 invention; 

FIG. 2 illustrates a diagram showing a current injecting device in FIG. 1 in 
accordance with a first preferred embodiment of the present invention; 

FIGS. 3 A and 3B illustrate diagrams each showing another embodiment of 
the current injecting device in FIG. 1 ; 
10 FIGS. 4 A and 4B illustrate diagrams each showing a current flow in the 

measuring object; and 

FIG. 5 illustrates a diagram showing another form of current flows in the 
measuring object. 

15 Best Mode for Carrying Out the Invention 

Reference will now be made in detail to the preferred embodiments of the 
present invention, examples of which are illustrated in the accompanying drawings. 
In describing the embodiments, same parts will be given the same names and 
reference symbols, and repetitive description of which will be omitted. FIG. 1 

20 illustrates a diagram of a system for visualizing conductivity and current density 
distribution in accordance with a preferred embodiment of the present invention. 

Referring to FIG. 1, the system for visualizing conductivity and current 
density distributions includes current injecting devices 100 for injecting currents in 
directions different from one another into a measuring object, such as a human body 
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or a substance, an MRI scanner 200 for measuring a magnetic flux density due to 
each of the currents, an operating part 300 for calculating a conductivity (or 
resistivity) distribution and a current density distribution inside of the measuring 
object from the currents and the magnetic flux densities induced thereby displaying 
5 means 400 for displaying the inside of the measuring object as an image of the 
conductivity and current density distributions, and a current source 500 for supplying 
the currents to the current injecting device 100. 

The current injecting device 100 includes an electrode 21 formed of a non- 
magnetic metal, such as copper, or a material having very high conductivity, and a 

10 wire for supplying the current to the electrode 21. Though direct attachment of the 
electrode 21 to a surface of the measuring object is possible, it is preferable that the 
electrode 21 is spaced apart from the surface of the measuring object for preventing 
distortion of the magnetic flux density at the surface of the measuring object. 
Therefore, as shown in FIG. 2, one embodiment of the current injecting device 100 

15 of the present invention further includes an insulating container 25 having the 
electrode 21 attached thereto for putting the electrode 21 on a place spaced apart 
from the measuring object. The insulating container 25 has one surface having the 
electrode 21 attached thereto, and the other surface to be brought into contact with 
the measuring object in an opened state. The insulating container 25 is stuffed with 

20 sponge having an electrolyte gel or electrolyte solution absorbed therein. The sponge 
having an electrolyte gel or electrolyte solution absorbed therein makes uniform flow 
of the current from the electrode 21 in the insulating container 25. Thus, by putting 
the electrode on a place spaced apart from the surface of the measuring object by 
using the insulating container 25, and supplying the current to the measuring object 
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uniformly by using the sponge having the electrolyte gel or electrolyte solution 
absorbed therein, distortion of an MRI image at the surface of the measuring object 
in contact with the electrode 21 can be reduced. 

Referring to FIGS, 3A and 3B, another embodiment of the current injecting 
5 device 100 of the present invention includes an electrode 31, a wire 33 connected to 
the electrode 31, an semispherical insulating container 35 having the electrode 31 
attached thereto. The insulating container 35 has one side having the electrode 31 
attached thereto, and the other surface to be brought into contact with the measuring 
object 'S' in an opened state. The insulating container 35 is stuffed with sponge 37 

10 having the electrolyte gel or electrolyte solution absorbed therein. The sponge 37 
may be replaced by only the electrolyte gel. 

The current injecting devices 100 are attached to a circumference of the 
measuring object 'S', and the current is made to flow only to one pair of the current 
injecting devices 100 at a time, such that the current flows through the inside of the 

15 measuring object 'S' with the one pair of the current injecting devices 100. If it is 
assumed that a number of the current injecting devices 100 is 'E\ a number C N' of 
cases of the pairs of the current injecting devices 100 for injecting the current into the 
measuring object 'S ? is E(E-l)/2. This is equal to a number of current paths in the 
measuring object C S\ A current between the pair of current injecting devices 100 is 

20 represented with I j (j = 1, 2, — , N). The current I J forms a current density J j (J j x, J J y, 
J J z) in the measuring object 'S\ 

The MRI scanner 200 measures a magnetic flux density B J in the measuring 
object 6 S' due to the current f, particularly, a z-directional magnetic flux density B J z 
parallel to a main magnetic field of the MRI scanner 200. That is, the MRI scanner 
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200 measures magnetic flux densities B ! z ? B 2 z ? B 3 z, — , B N z due to currents I 1 , 1 2 , 1 3 , 
— , I N injected into the measuring object 'S\ respectively. 

The operating part 300 selects one of the pairs of the current injecting 
devices 100 in succession, and controls the one of the pair of the current injecting 
5 devices to inject the current into the measuring object C S\ The operating part 300 
also applies the magnetic flux density B J z measured by the MRI scanner 200 to a 
preset algorithm, and calculates a conductivity a and a current density J j inside of the 
measuring object C S\ Since the conductivity a is an inversion of the resistivity p, the 
resistivity p can be represented with l/o\ Then, the displaying means 400 visualizes 
10 the inside of the measuring object C S 5 according to the conductivity a and the current 
density J j . 

A method for visualizing conductivity and current density distributions in 
accordance with a preferred embodiment of the present invention, and an algorithm 
applied thereto will be described, in detail. 

15 

FIRST EMBODIMENT 

Referring to FIG. 1, a plurality of current injecting devices 100 are attached 
to a circumference of a measuring object *S\ i.e., a human body, and the measuring 
object C S 5 is placed in an MRI scanner 200. In this instance, as described before, an 
20 insulating container 25 of the current injecting device 100 is stuffed with sponge 
having electrolyte gel or electrolyte solution absorbed therein. 

Pairs of the current injecting devices 100 are selected one by one in 
succession, and a current V (j = 1, 2, — , N) is supplied to an inside of the measuring 
object 'S' through the pair of the current injecting device 100 selected in succession. 
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For an example, as shown in FIGS. 4A and 4B, in a case four current injecting 
devices 100a, 100b, 100c, and lOOd are attached to a circumference of the measuring 
object 'S\ the currents are supplied to the inside of the measuring object 'S 5 in two 
directions in succession. At first, a current I 1 is supplied from the current injecting 
5 device 100a to the current injecting device 100c, and then, a current I 2 is supplied 
from the current injecting device 100b to the current injecting device lOOd. In these 
instances, the MRI scanner 200 measures the magnetic flux densities B l z and B 2 z 
due to the currents I 1 and I 2 injected into the measuring object C S' respectively. The 
MRI scanner 200 measures c z 5 direction components of the magnetic flux densities, 
10 only. 

Then, the operating part 300 applies the magnetic flux densities B*z and B 2 z 
measured by the MRI scanner 200 and voltages V J |as due to the currents I 1 and I 2 at a 
surface of the measuring object 6 S' to a preset algorithm, and calculates the 
conductivity a and the current density J j inside of the measuring object e S\ 
15 Equations to be used for the algorithm in the first embodiment will be 

described, where, 

S : a measuring object, 

dS : a surface of the measuring object, 

g 1 : a component of a current density of an injected current at a surface of the 
20 measuring object, 

n : a unit normal vector at the surface of the measuring object, 
V J : a voltage of the injected current in the measuring object. 
A relation between the current density J J = (J 3 x, J J y, J J z) and the magnetic 
flux density B j = (B J x, B j y, B J z) can be expressed as the following equation (1) 
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according to the Ampere's law J y = — VxB J . 

Mo 

- axJ J y+5yJ j x=— (5 2 x-i-a 2 y+a 2 z)B j z= — AB j z 

Mo Mo 

— (1) 

The equation (1) expresses a relation between the current density J J and the 
5 magnetic flux density B j z inside of the measuring object. In this instance, of the 
magnetic flux density B J , only a V direction magnetic flux component B J z is taken 
into account. 

J j -n = -cyVV j .it = g i 

— (2) 

10 The equation (2) expresses a relation between the current density and the 

surface voltage at a surface of the measuring object 'S\ 

The current density J J can be expressed as the following equation (3). 

ji = ^VV* 

— 0) 

15 Once the conductivity and the voltage inside of the measuring object are 

calculated, the current density of the injected current can be calculated by the 
equation (3). If equation (3) is substituted into a left side of the equation (1), the 
following equation (4) is obtained. 

SxadyV j - dyaSxV j = — AB J z 

Mo 

20 — (4) 

The equation (4) expresses a relation between the voltage and the current 
density inside of the measuring object C S\ As shown in FIGS. 4A and 4B, when two 
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15 



currents I 1 and I 2 are injected into the measuring object fi S' in different directions, the 
relation between the voltage and the current density inside of the measuring object 
'S 5 can be expressed as follows. 



'dyV 1 -dxV 
K dyV 2 -dxV 2 ) 



( dxa 









1 


AB 2 2 





- (5) 

In a case the two currents I 1 and I 2 are supplied to the measuring object 'S' 
through two pairs of the current injecting devices 100, if the current densities due to 
the two currents I 1 and I 2 satisfy the following equation (6), the conductivity a can 
be obtained. 



■j\ ji 



X 

■J 2 J 2 



*0 



— (6) 

Since the currents I 1 and I 2 of different directions are injected into the 
measuring object C S', the current densities due to the two currents I 1 and I 2 have 
different directions in almost all regions of the inside of the measuring object C S\ 
Therefore, the equation (6) is satisfied in all positions of the measuring objects C S\ 

If an inverse matrix of the equation (5) is taken, the following equations (7) 
can be obtained. 

9jcV 2 AB* -axV'AB* 
dxa = — lA — n z — > and 



dyer 



ju^dyV'dxV 2 -dxV'dyV 2 ) 

= dyV 2 AB 1 z -dxV 1 AB 2 z 
" Mfy^dxV 2 -dxV'dyV 2 ) 



20 



(7) 



WO 2004/062464 



13 



PCT/KR2003/002825 



Since the current density satisfies an equation d x 3 J x + d y J J y + 8 2 J J Z = 0 , the 
following equations (8) can be derived from the equations (7). 

V-((cr(x 0 ,y 0 ,z) + F[V\V 2 ,B z ,B z ])VV 2 ) = 0, 
5 (cr(x 0 ,y 0 ,z) + F[V\V 2 ,B z ,B z ])VV 1 -n as =-g\ 

(<t(x 0 ,y 0 , z) + F[V l ,V 2 ,B Z ,B 2 Z ]) V V 2 • n dS =-g 2 , 

JV 1= 0, and jV 2 =0 

-- (8) 

In the equations (8), the first and second equations are ones for the inside of 
10 the measuring object 'S\ and rest of the equations (8) are ones for the surface of the 
measuring object 'S\ F[V l ,V 2 ,B Z ,B Z ] can be expressed as follows. 

v f dyV 2 AB l z -dyV l &B 2 z 



l { uJdvY l dxV 2 -dxV'dvY 2 Y { 9 9 ) 



" MoidyY'dxV 2 -dxV l dyY 2 y 
The algorithm of the first embodiment calculates the conductivity a and the 
current density J j of an inside of the measuring object 'S' by using the foregoing 
15 equations. A process for calculating the conductivity a and the current density J J of 
an inside of the measuring object 6 S' by applying the magnetic flux densities B l z and 
B 2 z, and the voltages V j |as due to the currents I 1 and I 2 at a surface of the measuring 
object C S' to a preset algorithm, in detail. 

An initial conductivity Go is set to unity or any desired value, and by 
20 substituting the initial conductivity a 0 and the current density component gj at the 
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surface of the measuring object 'S' into the following equations (9), initial voltages 
Vo 1 and Vo 2 induced at the measuring object *S J by the injected currents I 1 and I 2 can 
be obtained. The initial voltages Vo 1 and V 0 2 are calculated by using the FEM or 
FDM in solving a partial differential equation, which are voltages for the initial 
5 conductivity ao set as desired. 

a Q (x, y, zyvVo -n\ ds = -g 1 , and 
J^-O.and 

10 c7 0 (x,j ? z)Vr 0 2 >n\ ds =~g 2 , and 

^ =0 

— (9) 

The initial voltages Vo 1 and V 0 2 calculated by using the equations (9), and 
the measured magnetic flux densities B l z and B 2 Z are substituted into the following 
15 equations (10), to calculate voltages Vi 1 and Vi 2 . 

V-((o-(x 0 ,y 0 , Z ) + F[F 0 , ,F 0 2 3 ^,^])VF 1 , ) = 0 5 

V • ((a(x 0 ,y 0 ,z) + F[V 0 \V 0 \B\,BlWV?) = 0 , 

(^ 0 ^ 0) z) + J P[F 0 1 ,r o 2 ,^, J B|])VF 1 , -n| 5S =-g 1 , 

{a(x, i y 0 ,z) + nV 0 \V 0 \B] z ,BlWV?-n\ ds =- g \ 

20 J>/=0, and [V? =0 

- (10) 
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Herein, a(x 0 ,y 0 ,z) + F\V* >B l z > B l ] > a value of the second conductivity 
&u can be obtained by using the initial voltages Vo 1 and V 0 2 , and the measured 
magnetic flux densities B* z and B 2 Z . The second conductivity ai is closer to a true 
conductivity than the initial conductivity a 0 > because the calculated conductivity cti 
5 has taken true magnetic flux densities B l z and B 2 Z of the inside of the measuring 
object 'S' into account. 

Then, the calculated voltages Vi 1 and Vi 2 , and the measured magnetic flux 
densities B* z and B 2 Z are substituted into the following equations (11), to calculate 

voltages V2 1 and V 2 2 . Herein, cr(x 0 ,y 0 ,z) + F[V l \V 1 2 ,B l z ,B%], a value of the third 
10 conductivity a 2 , can be obtained by using the calculated voltages Vi 1 and Vi 2 , and 

the measured magnetic flux densities B^ and B 2 Z . The third conductivity a 2 is closer 

to the true conductivity than the conductivities c?o and G\. 

Hereafter, by using the voltages V2 1 and V2 2 , a fourth conductivity CF3 and the 

next stage voltages may be calculated. 
15 V.((cr( Wo ,z)^^ 

V-((cr(x^*) + ^ 
(cr(Wo,*)+^^ 

J^-O.and j> 2 2 =0 

20 — (11) 

As described, for obtaining a conductivity close to a true value, values up to 
V^n+i and V 2 m+ i are obtained from the initial voltages Vo 1 and V 0 2 by using an 
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iterative algorithm. The calculation for obtaining the voltage values are repeated until 
the following equation (12) is satisfied. 
I V l m+ i- V^Kocuand 

I V 2 m+ i- V 2 m |<a 2 , (m = 0, 1, 2, 3, — -.) 

5 - (12) 

In other words, differences between the voltages at an (m)th step and the 
voltages at an (nrfl)th step are calculated, and it is determined whether the 
differences are smaller than oci and ot2 5 respectively. If the differences are greater than 
preset values respectively, the voltage calculation process is repeated. Opposite to 
10 this, if the calculated differences are smaller than the preset values respectively, the 
values of the V x m and V 2 m are substituted into the following equation (13), to 
calculate a conductivity a m +i. 

a = Umcr m =Um(a(x 0 ,y 0 ,z) + F[VlV*X,B,]) 

m~>oo m-»oo 

- (13) 

15 Even though 'm* approaches to infinitive in the equation (13), it is required 

that the 'm' is set to an appropriate value. The 'm 3 in the equation (13) increases until 
the difference of calculated conductivity satisfies a condition a m+ i-a m <8 . The 8 
is a set value. 

Then, a voltage corresponding to the calculated conductivity a m+1 is 
20 calculated by using the following equation (14). 
V.(aVF 1 ) = 0, 

Jf'-O. and 
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V-(oVF 2 ) = 0, 
oW 2 -n\ ds =-g 2 ,and 

Then, the voltages calculated by the equation (14) are substituted into the 
5 following equation (15), to calculate current densities. Then, by using the calculated 
conductivity a m+x or the current density, the inside of the measuring object can be 
visualized. 

J 1 =-oVV l ,and 

j 2 =-oVF 2 

10 (15) 

SECOND EMBODIMENT 

A plurality of current injecting devices 100 are attached to a circumference 
of a measuring object C S\ i.e., a human body, and the measuring object C S 5 is placed 

15 in an MRI scanner 200. In this instance, as described before, an insulating container 
25 of the current injecting device 100 is stuffed with sponge having electrolyte gel or 
electrolyte solution absorbed therein. 

Pairs of the current injecting devices 100 are selected one by one in 
succession, and a current f (j = 1, 2, — , N) is supplied to an inside of the measuring 

20 object 'S' through the pair of the current injecting device 100 selected in succession. 
For an example, as shown in FIG. 5, in a case four current injecting devices 100a, 
100b, 100c, and lOOd are attached to the circumference of the measuring object C S% 
the current injecting devices can be paired in six cases for injecting the currents, such 
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that the currents flow in six paths, and in six directions respectively under the control 
of the operating part 300. The currents of the six paths are a current I 1 between the 
current injecting devices 100a and 100c, a current I 2 between the current injecting 
devices 100b and lOOd, a current I 3 between the current injecting devices 100a and 
5 100b, a current I 4 between the current injecting devices 100b and 100c, a current I 5 
between the current injecting devices 100c and lOOd, and a current I 6 between the 
current injecting devices lOOd and 100a. In these instances, the MRI scanner 200 
measures the magnetic flux densities B*z, B 2 z, — , B 6 z due to the currents I 1 , 1 2 , — , 
I 6 injected into the measuring object 'S ? respectively. The MRI scanner 200 measures 

10 c z' direction components of the magnetic flux densities, only. 

Then, the operating part 300 applies the magnetic flux densities B*z, B 2 z, — , 
B 6 z measured by the MRI scanner 200 and voltages V j | 5S due to the currents I 1 , 1 2 , — 
-, I 6 at a surface of the measuring object 'S' to a preset algorithm, and calculates the 
conductivity a and the current density J j inside of the measuring object *S\ 

15 Equations to be used for the algorithm in the second embodiment will be 

described. 

It is assumed that the conductivity a inside of the measuring object C S' is 
isotropic within a range 0<a<oo. Since the conductivity a is not relevant to the 
current substantially, the voltages V J can be obtained from the following equation 
20 (16). 

V.(oVF y ) = 0,and 

-oVJ^.h = -V 0 = 1,2,3,-— ,N) 

~ (16) 

The current density g 1 is zero (0) at the surface of the measuring object C S' 
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10 



15 



not in contact with the selected pair of the current injecting devices 100. If the 
conductivity a, the injected current V, and positions of the current injecting devices 
100 are known, the voltage V j can be calculated by a calculating method, such as 
FEM. 

Moreover, an equation V 2 B = -// 0 VVx Vcr can be derived from the 

Ampere's law. If a z component of the equation is expressed with respect to a 
position of the inside of the measuring object 6 S\ the following equation (17) is 
obtained. 



— V 2 B{ = 



Mo 



dV J 



dV J 
dx 



da da 

dx' dy > 



(17) 



Where, uo is a magnetic permeability in the free space. The equation (17) 
can be expressed in a matrix form as the following equation (18). 



2 t>N 



dV l 


3V 1 


dy 


dx 


dr N 


dV N 


dy 


dx 



da 

dx 
da 

dy 



(18) 



If 



Mo 



V 2 ^ 



is set to be 'b', 



dv 1 


dV 1 


dy 


dx 


dV N 


dV N 


dy 


dx 



is set to be C U\ and 



do 

dx 
da 



is set to be 's', the equation (18) can be expressed as b=Us. 



WO 2004/062464 



20 



PCT/KR2003/002825 



If two currents I 1 and I 2 are supplied to the measuring object C S' through two 
pairs of the current injecting devices 100, if two voltages VI and V2 due to the two 
currents I 1 and I 2 satisfies the following equation (19), above matrix equation can be 
solved. 

dv 1 dv 2 av l av 2 A 

5 1 ^ 0 

dy dx dx dy 

— (19) 

Since the two current densities J 1 and J 2 due to the two currents I 1 and I 2 
have directions different from each other in almost all regions of the inside of the 
measuring object 'S\ the equation (19) can be satisfied in all positions of the 
1 0 measuring obj ect ' S 5 . 

The equation (18) can be expressed as the following equation (20) by using 
weighted regularized least square method. 

s -(U r U + AI)- i: U T b 

- (20) 

15 Where, U T denotes a transpose of 'U\ and U = WU, b=Wb, W=diag 

(wi, w 2 , — , Wn). W denotes an NxN matrix, I denotes a 2x2 unit matrix. X denotes 
a positive regularization parameter. The weighting factor Wj is set as the following 
equation (21). 

SNR f 

w j=~ " 

20 — (21) 

The SNR denotes a signal-to-noise ratio of the magnetic flux density B J Z 
measured by the MRI scanner 200. By adjusting the weighting factor Wj according to 
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the noise of the magnetic flux density B j z , a more accurate 's' value can be obtained. 

If the equation (20) is calculated for all positions of inside of the measuring 

IT 



objects 'S',a distribution of s = 
known. 



da da 

_ dx dy _ 



inside of the measuring object 6 S' can be 



For displaying the three dimensional (or solid) measuring object 6 S ? on a two 
dimensional screen, it is required to fix one of the three dimensions. Therefore, when 
it is intended to display a two dimensional image for convenience, it is assumed that 
z=0. It is also assumed that the conductivity at a fixed position r 0 =(xo ? yo, 0) at the 
surface 8S is unity. There are two methods for calculating the conductivity from 



10 V<7 = 



r da da 



, . One is the line integral method, and the other is the layer 

^ dx dy J 



potential method. 

The line integral method is the same with the method suggested in the first 
embodiment, in which the following equation (22) can be used. 

da , -y 7 & da t 



a(x,y) = a(x Q ,y Q ) + f o ^( x >^o>&+ j^^T^'J^ 



dx dy 
15 - (22) 

Where, a(x 0 ,y Q ) is set to unity. 

The layer potential method uses the following equation (23). 

a(r)={v 2 0(r~r'Mr')rfr' 

= - [ V r .<X>(r - r») • Vcr(r')dr'+ £ n r , • V r ,<D(r - r')a(r')dl r , 

- (23) Where, 0(r-r') = — loglr-r'L and V r .O(r - r ; ) = ~— . r ~ r ' . 

20 The following equation (24) can be derived with respect to the surface dS of 
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the measuring object 'S\ 

Urn j^ r , .V r ,®(r-tn r -rW}dl, = [«,'V r ,0(r-r'MrH - 

- (24) 

Moreover, the following equation (25) can be derived from equation (23). 

£^-±r^„ w= l((qlM2 lir . 



2;r ■° sr |r — r'| 2^ •» |r — r 

- (25) 

Where, eras denotes the conductivity limited to dS of the measuring object 
C S\ By using the equation (25), a surface conductivity eras, i.e., the conductivity at a 
surface dS of the measuring object 'S' can be obtained. By substituting the surface 
10 conductivity das into equation (23), the conductivity or of the inside of the measuring 
object 'S' can be obtained. By substituting the surface conductivity eras into equation 
(23), the following equation (26) can be obtained. 

^ = -|V r .O(r^/-0-V^)^+£^ -V,^^')^(r'H, — 

- (26) 

15 The equations (18), (22), and (26) can be repeated with respect to all regions 

of the measuring object 'S\ 

If the magnetic flux density B j z measured by the MRI scanner 200 is 
substituted into the equation (18), the 's', i.e., Va can be derived. However, if the 
voltage V j inside of the measuring object 'S* is not known, a true conductivity a 

20 inside of the measuring object 'S' can not be obtained. Therefore, an initial 
conductivity a 0 is assumed, to calculate the voltage V J inside of the measuring object 
C S\ and conductivity a m is calculated repeatedly until the conductivity a m 



WO 2004/062464 



23 



PCT/KR2003/002825 



approaches close to a true conductivity. In order to solve a problem in which the 
conductivity cr m differs from the true conductivity by a multiplication of a constant, 
the surface voltages V*)qs due to the injected voltage I J is used. The surface voltages 
V j |as are voltages at surfaces between the current injecting devices 100 that are not 
5 selected and the measuring object £ S' ? when the current is supplied to the measuring 
object *S* through the pair of the current injecting devices 100. The surface voltages 
V j |as are measured with the current injecting devices 100 that are not selected. 

A process for calculating the conductivity and the current density at the 
operating part 300 by applying the measured conductivity B j z and the measured 
10 surface voltages V j |as to the foregoing equations will be described in detail. 

The initial conductivity a m (m=0) is set, and the initial conductivity cr m 
(m=0) and the current density component g* at the surface of the measuring object 'S ? 
are substituted into the following equation (27), to calculate voltages V J m +i and 
V j m +i|as at the inside and the surface of the measuring object 'S\ Those are voltages 
15 for the initial conductivity set as desired. The current density component g 1 can be 
known from the current I J (j = 1, 2, — 9 N). 
V<a m VV j m+ 0 = 0, and 

-a m VVUi|as-n = ^ 0=1, 2, — , N) 

- (27) 

20 After the voltage V J m+i inside of the measuring object 'S' for the initial 

conductivity is obtained, the voltage V J m +i and the measured magnetic flux density 
B j 2 are substituted into the equation (18), and the equation (22), or equation (26) is 
used, to calculated the conductivity a m +i. The calculated conductivity a m +i is closer 
to a true conductivity than the conductivity a m , because the calculated conductivity 
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cJm-H has taken an true magnetic flux density B J Z inside of the measuring object 'S 5 
into account. 

However, even if the calculated conductivity a m +i has taken the true 
magnetic flux density B J Z into account, the calculated conductivity a m+ i can differ 
5 from the true conductivity. That is, though the calculated conductivity cr m +i has a 
value varied with the magnetic flux density B J Z , the calculated conductivity a m+ i may 
be greater or smaller than the true conductivity by a multiplication of a constant. For 
eliminating the difference caused by the multiplication of the constant of the 
calculated conductivity a m +i from the true conductivity, in the present invention, the 
10 measured surface voltage V J |as and the calculated surface voltage V m+ i j |as are used 
for scaling the conductivity o m +i. In other words, according to a ratio of the 
measured surface voltage V J |as to the calculated surface voltage Vm+i J '|as» the 
conductivity a m+ i is multiplied or divided by a constant value. 

Moreover, in the present invention, for obtaining the conductivity close to 
15 the true value, an iterative form of algorithm is used. That is, alike the following 
equation ||<x m+1 - a m || < s , a difference between a conductivity a m in a prior step and 

a conductivity a m +\ in the present step is calculated, and it is determined whether an 
absolute value of the difference is smaller than a preset value V or not. The equation 
(27) is calculated repeatedly by increasing the 'm' until the calculated difference 
20 becomes smaller than the preset value c s' . 

If the absolute value of the difference of the conductivity cr m in the prior step 
and the conductivity a m+ i in the present step is smaller than the preset value V, the 
conductivity a m+ i is taken as the true conductivity, by using the conductivity a m +i, 
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the measuring object *S* is visualized. 

Then, the conductivity a m+ i is substituted into the equation (16) to calculated 
a voltage for the conductivity a m +i, i.e., the true voltage V 3 inside of the measuring 
object *S' is calculated. Then, by using the conductivity a m+ i and the voltage V J , the 
5 current density J j is calculated, and by using the current density J J , the inside of the 
measuring object C S' is visualized. 

Industrial Applicability 

As has been described, the present invention permits to obtain a conductivity 
10 distribution image and a current density distribution image of an inside of a 

measuring object without rotating the measuring object, such as a human body, or a 

substance, within an MRI scanner. 

Moreover, by using the current injecting devices together with the present 

MRI scanner, the conductivity distribution image and the current density distribution 
15 image of the inside of the measuring object are made available, that have not been 

available only with the present MRI scanner, and by using which a high quality 

image can be reproduced. 

Furthermore, the conductivity distribution image and the current density 

distribution image can be applicable to the present medical equipment. According to 
20 this, a functional imaging can be provided for organs each of which biological tissue 

conductivity varies with physiologic action, to secure a new medical examination 

technology. For an example, an image required for measuring a volume of a lung 

varies with respiration (a measurement by using variation of a resistivity of the lung 

in inhale and exhale of air), a measurement of an cardiac output of a heart, a 
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measurement of a stomach, a measurement of cystic residual urine, a measurement of 
a healed state of a bone fracture, or a measurement of a brain function. 

Since quantitative electric properties of a human tissue can be known, 
optimization of diagnosis equipments, and treatment equipments using electro- 
5 magnetic energy, such as various electric/magnetic stimulators, defibrillators, cardiac 
pacemakers, and the like, is made possible. 

It will be apparent to those skilled in the art that various modifications and 
variations can be made in the present invention without departing from the spirit or 
scope of the invention. Thus, it is intended that the present invention cover the 
10 modifications and variations of this invention provided they come within the scope 
of the appended claims and their equivalents. 
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What is Claimed is: 

1. A system for visualizing conductivity and current density distributions 
comprising: 

a plurality of current injecting devices for injecting currents into a measuring 

5 object; 

an MRI scanner for measuring one directional component of a magnetic flux 
density due to each of the currents injected into a measuring object; 

an operating part for controlling the current injecting devices so as to inject 
currents of different directions into the measuring object, and calculating a 
10 conductivity distribution and a current density distribution inside of the measuring 
object by using the one directional component of a magnetic flux density; and 

displaying means for visualizing the conductivity and current density 
distributions calculated by the operating part. 



15 2. The system as claimed in claim 1, wherein the current injecting device 

includes; 

an electrode, 

an insulating container with the electrode attached to one side, the insulating 
container having an electrolyte substance, and 
20 a wire for supplying the current to the electrode. 



3. The system as claimed in claim 1, wherein the operating part controls the 
current injecting devices such that one pair of the current injecting devices are 
selected in succession, and the selected pair of the current injecting devices supply 
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the current to the measuring object. 

4. The system as claimed in claim 1, wherein the operating part calculates an 
inside voltage and a surface voltage of the measuring object for a first conductivity, 

5 and calculates a second conductivity by using the inside voltage and the one 
directional component of the magnetic flux density. 

5. The system as claimed in claim 4, wherein the operating part multiplies or 
divides a constant to the second conductivity according to a ratio of a calculated 

10 surface voltage to a measured surface voltage. 

6. The system as claimed in claim 4, wherein, if an absolute value of a 
difference of the first conductivity and the second conductivity is greater than a 
preset value, the operating part calculates the inside voltage and the surface voltage 

15 for the second conductivity and calculates a third conductivity by using the inside 
voltage for the second conductivity and the one directional component of the 
magnetic flux density. 

7. The system as claimed in claim 4, wherein the operating part determines 
20 that the second conductivity is a true conductivity if the absolute value of the 

difference of the first conductivity and the second conductivity is smaller than the 
present value. 

8. The system as claimed in claim 4, wherein the operating part substitutes 
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the inside voltage and the one direction component of the magnetic flux density into 
the following equation, and subjects the equation to a line integral, to obtain the 
second conductivity. 



Mo 



V 2 B" 



dV 1 


dV 1 


dy 


dx 


dv N 


dV N 


dy 


dx 



da 

dx 
da 

ldy_ 



Where, j^o denotes a magnetic permeability of the free space. 



10 



9. The system as claimed in claim 4, wherein the operating part substitutes 
the inside voltage and the one directional component of the magnetic flux density 
into the following equation, and solves the equation by layer potential method, to 
obtain the second conductivity. 

dV 1 dV l 



Mo 



V 2 B X 




da 

dx 
da 



Where, \xq denotes a magnetic permeability of the free space. 



10. A method for visualizing conductivity and current density distributions 
15 comprising the steps of: 

(a) injecting currents of different directions into a measuring object through 
current injecting devices attached to a surface of the measuring object; 

(b) measuring one directional component of a magnetic flux density of due 
to each of the currents injected into an inside of the measuring object; 
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(c) calculating a conductivity and a current density of the inside of the 
measuring object by using the measured one directional component of the magnetic 
flux density; and 

(d) visualizing the conductivity and the current density. 

5 

11. The method as claimed in claim 10, wherein the step (a) includes the 

steps of; 

selecting a pair of the current injecting devices in succession, and 
injecting the currents into the measuring object in succession through the 
1 0 selected pair of the current inj ecting devices . 



12. The method as claimed in claim 10, wherein the step (c) includes the 

steps of; 

calculating the inside voltage and the surface voltage of the measuring object 
15 for a first conductivity, and 

calculating a second conductivity by using the calculated inside voltage and 
one directional components of the magnetic flux density. 



13. The method as claimed in claim 12, wherein the step (c) further includes 
20 the step of multiplying or dividing the second conductivity by a constant according to 
a ratio of the calculated surface voltage and a measured surface voltage of the 
measuring object. 



14. The method as claimed in claim 12, wherein the step (c) further includes 
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the steps of; 

calculating an absolute value of a difference of the first conductivity and the 
second conductivity, and 

comparing the absolute value of the difference to a preset value. 

5 

15. The method as claimed in claim 14, wherein the step (c) further includes 
the steps of; 

calculating the inside voltage and the surface voltage of the measuring object 
for the second conductivity if the absolute value of the difference of the first 
10 conductivity and the second conductivity is greater than a preset value, and 

calculating a third conductivity by using the inside voltage for the second 
conductivity and the measured one directional component of the magnetic flux 
density. 

15 16. The method as claimed in claim 14, wherein the step (c) further includes 

the steps of determining the second conductivity being a true conductivity, if the 
absolute value of the difference of the first conductivity and the second conductivity 
is smaller than the preset value. 

20 17. The method as claimed in claim 12, wherein the step of calculating a 

second conductivity by using the calculated inside voltage and the measured one 
directional components of the magnetic flux density includes the steps of; 

substituting the inside voltage and the one directional component of the 
magnetic flux density into the following equation, and 
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subjecting the equation to a line integral, to obtain the second conductivity. 

'da 



Mo 



2 tdN 



V 2 B 



dV 1 


dV 1 


dy 


dx 


dV N 


dv N 




dx 
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da 

dy 



Where, jliq denotes a magnetic permeability of the free space. 
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18. The method as claimed in claim 12, wherein the step of calculating a 
second conductivity by using the calculated inside voltage and the measured one 
directional components of the magnetic flux density includes the steps of; 

substituting the inside voltage and the one directional component of the 
magnetic flux density into the following equation, and 

solving the equation by layer potential method, to obtain the second 
conductivity. 



1 






dV 1 
dy 


dV 1 

dx 


Mo 


V 2 5f 




dV N 


dV N 






dy 


dx 



da 

dx 
da_ 

dy 



Where, |Uo denotes a magnetic permeability of the free space. 



15 19. The method as claimed in claim 10, wherein the step (c) includes the 

steps of; 

calculating an initial inside voltage for any preset initial conductivity, 
calculating a conductivity by using the calculated initial inside voltage and 
the measured one directional component of the magnetic flux density, and 
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calculating the current density of the inside of the measuring object by using 
the calculated conductivity and the inside voltage of the measuring object for the 
conductivity. 

5 20. The method as claimed in claim 19, wherein the step (c) further includes 

the step of comparing the inside voltage in a previous step and the inside voltage of 
the measuring object for the conductivity. 
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FIG. 2 
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FIG. 4B 
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